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ABSTRACT 


This thesis presents a computer program designed to 
estimate the first seven natural frequencies and their 
respective modes for ships and to estimate ship response 
in a Pierson-Moskowitz sea. The program was designed to 
function on a minimum of input. Hull offsets and deck edge 
height for each station are the major data required. The 
program computes the hydrodynamic added mass, the polar 
moment of inertia J, the second moment of area I, the 
effective shear stiffness KAG and the ship mass for each 
station. A Prohl sequence is used to determine hull 
natural frequencies and mode shape. A wave forcing function 
is presented which sums buoyancy, added mass and damping 
factors due to wave passage. The program can be used for 
all vessels, but program development was based on a single 
beam model so superstructure effects are not considered. 
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a. INTRODUCTION 
Vibration is a phenomenon specially appealing to a 
number of engineers due partly to a clearly defined, exten- 
Sive and intriguing body of underlying mathematical analysis. 
From the late 1800's when attention was first paid to the 
problem of ship hull vibration this body of analysis has been 
expanded by a large number of active researchers to provide 
the current comprehensive level of understanding of the subject. 
It is apparent that vibration of the ship hull is a 
source of concern for the naval architect (and shipowner). 
A passenger vessel with vibration amplitudes above a certain 
threshold would have few repeat passengers. Fatigue of 
structural components due to the cyclic stresses is a factor 
to be considered. The conventional method of determining the 
Strength of the hull girder does not explicitly take into 
Mecount the cyclic stress associated with hull vibration. 
It would seem necessary to determine what these stresses 
might be, how they vary with time and the effect of their 
addition to stresses already present through ship loading 
and the quasistatic effect of wave support. The sophisticated 
equipment on board warships experience problems in an environ- 
ment of excessive vibration. It is interesting to note that 
with continued refinement of such equipment as navigation 
instruments, radar, sonar, fire control computers, communica- 


tions gear and electronic warfare devices the requirement for 





vibration free foundations becomes more stringent. There 

are dangers to shipping from acoustic and pressure operated 
offensive weapons. Attention has been focused on noise 
emitted by hulls and hull vibration in general. The "acoustic 
Signatures" of vessels play an important part in tracking 
foreign vessels, both surface and subsurface, by passive means. 
It appears that the problem of vibration control is becoming 
more important. 

As stated, the body of analysis is large but it is 
apparent that the complexity of hull structures in both form, 
material and construction and the complexity of three-dimen- 
Sional hydrodynamics combine to limit the application of the 
analysis to reality. It is not surprising then that the pro- 
blems of hull vibration still remain and that work continues 
on new methods of mathematical treatment, on models and on 
full size ships where the adequacy of the technical progress 
can be shown. 

This thesis is concerned with one of the two broad 
Classes of ship hull vibration, that is synchronous or 
resonant vibration as opposed to local vibration. In this 
case the entire hull acts like a beam vibrating under some 
excitation from machinery, the propellor or waves. [In 
serious cases the motion can clearly be seen by sighting 
along the length of the ship. The importance of a resonant 


Vibration is dependent on the conditions of its occurence. 





T£ the engine at normal operating speed excites the hull then 
a serious problem exists. If, however, this resonance occurs 
at an engine speed that will be only infrequently maintained, 
then there is a problem of considerably less significance. 
Similarly, if a vessel is expected to operate in certain 
known sea states then it is necessary to avoid an excessive 
resonant condition by proper hull design. 

The other class of hull vibration is a local effect. 
Here specific parts of the ship or certain fittings such as 
masts, stacks, superstructure or a panel of plating are set 
into motion. Depending on the point of view local vibration 
can also include the undesired motion of items as small as 
ment fLixtures or stanchionchains. Local vibration can have 
the same effects as hull resonance, that is personal discom- 
fort and structural and equipment failures. This type of 
Vibration is not as readily predicted and analyzed as whole 
hull vibration and is treated generally after the fact by 
locating and removing the excitation or adding local 
stiffening, pillars and similar structures. Obviously in 
many cases the excitation cannot be eliminated and additional 
Structure is not feasible. Then the solution is the removal 
of the object or relaxation of the requirements. 

It is apparent that the hull of a ship will vibrate 
only if an external force is applied to it. When this force 


is from slamming, the resulting shudder can be picked up by 





accelerometers and can be shown to be transient in nature. 

A similar force can, be generated by rapidly releasing a 
haused anchor or rapidly braking a free falling one. The 
transient response and information which can be obtained or 
it will be discussed later. Other sources of excitation can 
produce continuous vibration which is somewhat more insidious 
due to its gradual deleterious action. Such sources include 
those whose origin is some out-of-balance periodic force in 
the diesel main engine, auxiliary machinery, shafting or 
propellers. 

It can therefore be reduced in any given instance by 
proper attention to balancing during construction. Diesel 
engine balancing is a well developed art which depends on the 
type of engine, number of cylinders, auxiliaries which may be 
run from the main engine and many other factors. Properly 
Operating turbine equipment runs at angular velocities well 
above the frequencies of concern here. 

The size of shafting and propellers makes perfect 
balancing impossible. In addition other once per revolution 
forces can be created if one blade of the propeller has a 
different shape than the other blades. Even with a perfect 
vropeller, the blades work in a mixed wake behind the hull and 
the force on each varies continuously throughout each revolu- 
tion. These varying forces are conveniently divided into two 


types according to the method of transmission to the hull. 





"Propeller" forces are transmitted through the shafting to 
struts and bearings. "Surface" forces are transmitted 
through the water to the stern. Both types have a frequency 
of propeller angular velocity times the number of blades and 
are referred to as blade-rate forces. Extensive work has been 
carried out and resultantly the magnitude of these forces can 
be calculated quite accurately. For example, the propeller 
forces can be determined from lifting surface theory assuming 
the hull induced velocity field in the propeller plane is 
known. With the requirement that the fluid particles should 
remain in contact with the propeller blades, it follows 

that the velocity components normal to the surface of the 
blades due to the hull induced wake should be equal to the 
corresponding components of the propeller induced fluid 
velocities. This requirement can be expressed as a three- 
dimensional integral equation. The numerical computation is 
formidable but computer solutions yield the pressure distri- 
bution and the blade lift. Integration yields the forces. 

It is interesting to note that due to the dynamic behavior 

of the shaft, bearing reactions can be an order of magnitude 
higher than the forces generated at the propeller. All of 
the forces can be kept within tolerable limits by careful 
manufacture and proper attention to stern and appendage 
design and particularly to blade tip clearance. Some 
hydrodynamic disturbing forces have to be accepted, however, 
as part of the price of mechanical propulsion. 
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The excitation due to wave passage is at present an 
active area of research. The circular frequency of encounter 
of waves, Wee is dependent on the wave spectrum of the sea 
surface and the speed of the vessel. Even for the swiftest 
of destroyers the amount of energy available from excitation 
of even the lowest natural frequency 1S very small, the 
preponderance of wave energy being at lower frequencies. 
Ships of more recent vintage, however, are longer with some 
lengths over 1,000 feet and natural frequencies for the two 
node mode of between 3 and 5 radians/sec. There is ample 
energy for excitation in this range and numerous studies have 
concluded that "springing" as two node vibration is called 
mona Significant addition to midship stress levels. [1] [12] 
m6} [30] 

When engines were first put on shipshull vibration 
problems arose. Shortly thereafter an effort was made to 
predict the frequency for the two node mode of vibration. 

As a result of very early work by Otto Schlich and by F.wW. 
Todd, empirical formulae were developed which closely 
correlated with observed frequencies. Schlich proposed a 
modification to the natural frequency formula for a uniform 
free-free beam which is 

(B51) 7 
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where BoL = 4,730, E is the Young's modulus, I is the section 
moment of inertia, m is the beam mass per unit length, and 
L is the length. This formula will be explained more fully 


later. Then 





oe) = N =a: 2) 


minere the units are tons for A, feet for L, et mor 2 and 


radians/sec for Wo Schlick gave the following values for N: 
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Another early effort which was twin steps in the 
directions of simplicity and accuracy by F.W. Todd was based 


on 


i 765) (3) 


where B is the beam, D is the depth, and C a coefficient 
dependent on L/D and the ship type. In the Todd formula, 
the virtual mass was used which was a function of the beam 


Moraratt ratio B/T. So 








: (4) 
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where 


i 


Palins 2a 3) (5) 


Values for g are listed in [32]. 


In the treatment which follows the natural frequencies 
of a vessel are determined from offset data and other infor- 
mation. The process iS an iterative one and can be expedited 
by an appropriate choice of a starting frequency. The choice 


@e this value is from a formula devised by L.C. Burrill. 


D = 80,0007 [——_—_4+___}/"" (6) 
ANI (1+B/2T) (1+r_) 
where 
2 3 2 
——— 3.5D” [3 (B/D) ~+9 (B/D) “+6 (B/D) +1.2] (7) 
S 


L* (3B/D+1) 


The Taylor shear correction term t . is dimensionless. 

The analysis of hull vibration has proceeded far beyond 
that which produced the preceding empirical formulations and 
this work attempts to explore the advances in considerations 


of damping, excitation, added mass and the use of the computer. 
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Zi, TECHNICAL TREATMENT 

A. Uniform Beam Vibration 

As previously stated, the mathematical base for the 
analysis of vibrating bodies is large and acceptance of 
certain broad assumptions allows complex problems to be 
idealized and working solutions obtained. As more realistic 
assumptions are made more intimate knowledge of the process 
is Obtained. The price,of course, is vastly more complex 
solution techniques. 

Consider a beam of length L. The vertical displacement 
at any time, t, is y(x,t) where x is the longitudinal axis 
of the beam. The total vertical force per unit length is 
£f(x,t). The system parameters are the mass per unit length 
m(x) and the flexural rigidity EI(x) where E is Young's 
modulus of elasticity and I(x) is the cross section area 
moment of inertia about the neutral axis of bending at any 
position x along the beam. Consider an element of length 
dx of the beam. V(x,t) and M(x,t) are the shearing force 
and bending moment. The simplifications to be included at 
this point are 1) shear deformation is small in relation to 
bending deformation and 2) rotary inertia effects are small. 

These simplifications are valid if the beam's length 
to depth ratio L/D is about 10 or greater and if deflections 
and slopes are small. Most ships underway fall into this 


category. [14] Another major assumption is that damping is 


eZ 





small enough to be neglected. This is sufficiently accurate 


moo Vibrating ships particularly in the free vibration 


x 


analysis which follows. [12] The small damping present 
oe 


serves to limit vibration amplitude rather than modify 
frequency or vibration pattern. The force equation of motion 


in the vertical direction is then 


2 
(v(x, t) +e) ax] - V(x,t)+£ (x,t) dx=m (x) dx2 Yt) 
x 542 
(1) 


The moment equation of motion, ignoring the inertia torque 


associated with the element rotation is 


oOM(x,t) 
Ox 


oV(x,t) 


ase axiax 


[M(x,t) + ax]-M(x,t)+[V(x,t)+ 


Pe (x,t) ax = = 0 (2) 


Here we assumed that the elemental length dx is small enough 
eemtnat £(x,t) 1s constant along the element. Equation (2) 


reduces to 


eb ( sa, tc) 


ae ov (x,t) = 0 (2) 


Using (3) in (1) and simplifying results in 


2 2 
0 ee EI A fice oye m(x) a yixet) (4) 
ox at 


i 





From basic strength of materials it is known that 


2 
M(x,t) = EI(x) oytet) (5) 
Ox 
(5) into (4) yields 
2 2 2 
~2 (er (x) 2 YS *) | 48 (x,t) z: m (x) 2 ee) (6) 
Ox Ox ot 


which is the differential equation for the flexural vibration 
of a beam. [25] This is a fourth order linear nonhomogenous 
partial differential equation. General solution theory re- 
quires a number of boundary conditions equal to the order of 
the equation. Since this equation must be integrated four 
times to get a solution, four arbitary constants are intro- 
duced, thus the requirement for four boundary conditions. 

In the case of interest the beam ends are free. Thus the 


Shear and moment at the ends equal zero. From equations (3) 


and (5) 

ace) 

M(x,t) = EI (x) cet = 0 Ga 
Ox 
5 yy (x,t) 

Wace) Sealer (x) See) = 6 (8) 

Ox 2 

ox 

ae —iae Ores iF 
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The next analytical step is consideration of free vibration. 


The equation of motion reduces to 


2 2 2 
Ul ‘) a 
a ort (x) .) (x sy = m(x) x7) (9) 
2 2 2 
Ox Ox dt 
with the boundary conditions (7) and (8). This forms a 


classic "boundary value" problem. It 1S opportune at this 
point to explore the possibility of synchronous motion. This 
is a motion in which the general shape of the beam displace- 
ment does not change with time. That is to say, every point 
of the beam executes the same motion, passing through the 
equilibrium position and the point of maximum excursion at 
the same time. This implies that the deformation g(x,t) is 
separable in space and time. Therefore a possible solution 


of this boundary value problem is of the form 
VRCne = (XxX) E(t) (10) 


where Y(x) represents the general beam configuration in 
vertical vibration and depends on x alone, and where F(t) 
indicates the type of motion the beam executes and depends 


@ot alone. Introducing (10) into (9) 


ant () 


Y(x)m(x) 2 (11) 
D4 


Ike 





where partial deriviatives have been replaced by total 

derivatives. In addition, the variableshave been separated 
so that the left hand side depends on x along and the right 
hand side on t alone. It is apparent, therefore, that both 


Sides must equal a constant. Let A’ be that constant, then 


ar (t) _ 


at? 


AF(t) = 0 (EZ, 


where, since Y(x), m(x) and I(x) are real, A is real. 


Let F(t) be of exponential form 


F(t) = Ae (13) 


mgen Substitution of (13) into (12) yields 


t t , 


Seon dividing by Ae~ (for Ae>* # 0) 


SO 


s = ty i (EA) ) 
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If A is a positive number then the two roots 


real and equal but opposite in Sign. 


(14) are both 


This suggests that 


there are two solutions for F(t), one increaSing and one 


decreasing exponentially with time. 


are inconsistent with an undamped conservative system so 


the possibility of A being positive can be discarded. 


Letting A = ane we have s = +1w and 


Expansion of this expression where 


cos wt + 1 Sin Wt 


1) 
Hl 


yields 


E(t) 1 1 


A.cos Wt + iA.sSin Wt + A.cos Wt - 


2 


(A, +A.) cos Wwe + 1(A,-A,) sin Wt 


iA,sin Wt 


This must reduce to a real expression and will if 


A, + A, = Ccosd 


i(A,-A,)= Csing 


Ev 


Both of these solutions 


G15) 





then 


Pie = CCOS tt = 9?) (16) 


Proceeding with the other part of (11) yields 


2 Z 
Ee = we mix) Y¥ (x) (17) 
D Z 
ax ax 
Oa only 
No general closed form solution of (17) exists. How- 


ever, for certain special cases solutions can be obtained and 
for illustrative purposes the following major simplification 
will be made, EI(x) and M(x) are made constant along the 
beam's length. This takes all ships with the possible ex- 
ception of some box shaped barges out of consideration but 


the information derived from such a step will provide some 


weight into the problem. (17) reduces under these assumptions 
15(8) 
4 2 
= BY(x) = 0 go = om (18) 
alike 


The general solution of (18) can readily be verified to be 


¥ (x) =A, Singx+A,cosgxtA Sinhgx+A ,coshgx des) 


2 - 
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Simplifying the boundary conditions (7) and (8) yield 


q-vin) ay 
aes 


a°x ass 


ape 


The second and third derivatives of Y(x) are 
hae = Bo. 2 2: Z 
Né (x) =A, 8 Sinfx-A,P cOSBxtA,8 sinhsx+A,p coshBx 
See fi - 3 Sor 3 ca 
Me (x) =A, 8 cosBxtA.,8 sin6x+A,8 cosh8xtA,6 sinhgx 


meex=O0 (20) and (21) yield 


At x=L substitution into (20) yields 


A. (Sinh&L-singL) +A 


1 (cosh@L-cos8L) =0 


2 


emer Substitution into (21) yields 


A. (coshgL-cos8L) +A 


l (Sinh8L+sin8L) =0 


2 


(21) 


(23) 


Parenthetically, these two equations in two unknowns can be 


irterms of A 


reduced to one equation leaving A, l 


aE 





-A. (cosh8L-cos8sL) 
ie 7 ia (25) 
Z sinh&Lt+sin8L 


Maserting (25) and (22) into (19) yields 


(cosh8L-cos8L) 


Sinhe@LtsingL (CoSbxtcoshséx) | (26) 


Y (x) =Af{sin8x+sinh8x- 
From linear algebra the determinant of the coefficients in 
(23) and (24) must equal zero for non-trivial values of Ay 
and A, to exist. That is 


2 


sinhfL-singL cosh8L-cos8éL | 
= at) 
cosh&L-cosgL SinhsL+sin8gL | 
mars is known as the "characteristic determinant" and reduces 


Ee, 


Eero sinhci sinemeicosheL-eoser)~ = 0 
After simplification this yields 
EQsnenCcosel = 1 (27) 


This is known as the "characteristic equation" and it is 
apparent that it can be solved for an infinite number of 
values of 8,that is, for oo where r=1,2,3... which are defined 
by 

20 





ge = we m/EL (28) 
From (28) an infinite set of frequencies, Wi», are obtained. 
These are natural frequencies. For each O equation (26) 
yields a prescribed vibration shape ee) known as a "natural 
mode." 

The general problem of determining the values of 8 in 
feo) £Or which nontrivial solutions Y(x) exist where certain 
boundary conditions are prescribed is called the "characteris- 
tic-value" or "eigen-value" problem with the parameters w 
being the eigenvalues and the functions ¥ (x) being eigen- 
functions. [19] [25] It is to be noted that the eigenfunctions 
are precise as to configuration but not to absolute value. 
Since (18) is homogeneous and ¥ (x) is a solution, then so 
is ay (x) where ® iS a constant multiplier. A unique property 
of these eigenfunctions forms the basis for this entire train 
of analytical thought. For the general case of beam vibration 
as expressed in (6) this property known as "orthogonality" 
1s expressed by 

L 

f m(x) Y (x) Ye (x) = 0 r#s (29) 

O 
This orthogonolity is with respect to the mass function which 


serves as a weighting term. 


ok 





This orthogonality implies that any possible configura- 
tion of the beam can be represented as a linear combination 


of the natural modes of the beam. That is 


y(x) = 2 C. Y (x) (30) 


where the C.. have been replaced by the time functions n(t) 
which are known as the "natural coordinates." 

As previously stated the eigenfunctions are not deter- 
Mined in regards to absolute value. A convenient normaliza- 


tion scheme is such that the eigenfunctions satisfy 

L 

d m (x) XY (x) Y¥, (x) = 6 Ce, ae (2) 
where oe PomenemGiracnGaelta Lunction 


0 r#s 


S sae z 
1 r=s 


In the uniform beam example insertion of one eigenfunction 
and eigenvalue into (18) yields 
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31> ar een Pome Y. Gx) (35) 


Multiplying both sides by Y, (x) and integrating over the 


length gives 


IE aty (x) L 3 
ower (x) EL = = {eae 3h (Bye (x) 
os dx O 


7 a ee i) 


Now substituting the assumed solution (31) into the beam 


equation (9) and assuming a uniform beam yields 


fore) co a*y (x) 

DS reo ame (x) + 2 q(t) EI———— = 0 (3:5) 
ite te ic 4 

r= ta— pb ax 


Multiplying through by Y 5 (x) and integrating over the length 
and considering (32) and (34) yields the independent set of 
ordinary differential equations 
h(t) + wr n(t) = 0 r=1,2,3 (36) 
ic gman a aS 
The solution to each of these equations can be shown to be 


Nn. = Cc. cos (wt ~ >) OS) 


where C.. and o. are determined from initial conditions of the 


Beam, y(x,o) and y(x,o). 
23 





Using (37) and (31) and setting t=0 yields 


y (X,0) C.. cos? Y (8) (38) 


| 
ides 


| 
4 


ix CO) tC _w sine, ¥. (x) 


=1 


Mottiplying by mY . (x) and integrating over the length yields 


L 
C.coso, = ; mY,(x) y(x,0) dx 
(2) 
Loin : 
Cc -sing, - ; we Y-(x) y(x,0) dx 


which can be solved for C.. and o-- 

The free vibration case for the uniform free-free 
beam has now been completely solved. The natural modes 
wee) and the natural frequencies w have been computed. 
It is the orthogonality of these modes and the fact that 
they span the vector space in which they are contained that 
makes the expansion theorems (30) and (31) correct. The 
natural coordinates n- (t) can be determined from initial 
@enditions. 

The forced vibration case is slightly more complex. 
Here substitution of the assumed solution (31) into the 


beam equation (6) yields 
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a 


ne mage nie eee (=) BE ——— = £ (x,t) (40) 

r=] 7 r=1 dx 
By multiplying by ¥ , (x) and integrating over the length and 
considering the normalization integrals (32) and (34) the 
following independent set of ordinary differential equations 


are obtained 


figs) © ues) = Tae) (41) 
where 
L 
F(t) = , Ys) 16 (ron ec laaeb. (42) 


The solution to (41) generally involves the use of the 


Senvolution integral 


c 
i) (ae) = 4 erect) Sin Wlt—-t) dT 
IG U) ele: pe 
ig ano) 
+ n- (0) cos wt - n- (9) Sin wt (a3) 


where n,(0) and n_(0) are as determined by (39). 

This completes the classical analysis of the forced 
vibration of the free-free uniform beam. It is not difficult 
Bemsee that significant complications arise when the desire 


is present for a closer representation of reality. [In (7) 
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inclusion of the beam theory damping term results in 


Rented) ja 
M(x,t) = EI(x) — + B(x) ene ann (7a) 
Ox ox oat 


where 8@(x) represents viscous structural damping. More on 


this in Section E. The equation of motion (6) then becomes 


Z Z 2 Z 3 
m (xX) oye t) + 2(EI (x) 22 *)) fe 2 (B(x) 3 pe 5) 
dt ote ax Ox ax Ot 
= £ (x,t) (6a) 


A free vibration analysis would proceed with the solution of 
(9) as before but no assumption as to uniformity would be 


made. The resulting general solution 


W(x) = AF(x) + BG(x) + CA(x) + DJ(x) (44) 


can be obtained for sufficiently well behaved EI(x) and m(x). 


Boundary conditions (7) and (8) require that 


Ao je teaG (Ol) CH" (0) + DI" (0) = 0 
at )ee ie BG (i) + CH (LL) + DJ" (L) = 0 
Pee CO mec eGo) + CH"(0) + DI™(@) = 0 
Ae ie teeGih) + CH(L) + DI'(L) = 0 
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This leads to the characteristic determinant 


F" (0) econ H" (0) 10) 
F{" (ty G" (in) H" tr) x" (i) 
F™(0) GLO } H'" (0) 0) 


Beer (iy) Cpt ils) Jel (IR) wien) 


This is an equation in w and its roots are the system 
natural frequencies. For each root W there correspond 
specific values of A, B, C, and D and therefore a specific 
eigenfunction ve: tae snature OF Lhe four solutions F, G, H, 
and J may make the solution of the determinant quite complex. 

This is indeed an understatement. The general solution 
(44) will be a combination of transcendental functions and 
infinite series. Solution of the determinant would be very 
time consuming and costly. Additionally, inherent in the 
previous work has been the assumption that shear defection 
and rotatory inertia are small factors. This is reasonable 
for the lowest vibration modes but increasingly becomes a 
source of error as the mode number increases, say for three 
or more nodes. The shear defection, for example, is about 
15% of the bending deflection for the lowest vibration mode, 
miat Of two nodes for L/D = 15. [14] For vibration of about 
four nodes shear deflection becomes greater than bending 


deflection. 
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Bere sone seCouchece 

A method which allows account to be taken of shear and 
rotary inertia and can be solved for the eigenvalues and 
functions is the technique devised by M.A. Prohl in 1945. 
It was originally produced for calculating the angular 
velocities at which flexible rotors would resonate. The 
actual rotor was simulated by a "lumped mass" system in 
which the rotor mass was concentrated at a number of stations. 
Station spacing was a matter of convenience though obviously 
the greater the number of stations the greater the accuracy. 
The stations were connected by weightless rods whose bending 
stiffness and shear characteristics were identical to those 
of the rotor section. This effective discretization allows 
problem solution by numerical techniques. The procedure is 
directly applicable to the analysis of ship hull vibration 


as will be shown. Let 


a = station mass plus added mass 
M; = bending moment 
Vs = shear force 


g; = vertical deflection 


Os = slope 
be = length between adjacent stations 
w = angular velocity 


KGA, = shear stiffness which is constant along L. 
EI. = bending stiffness which is constant along L. 
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Uo =sass moment Gf inertia 


subscript denoting consecutive stations. 


- 
I 


stayel 


The equatiors from strength of materials which govern the 
motion of the stations and which include shear and rotary 


inertia terms are 


9 ee pence Vge1bi+12 (1) 
itl aL Eley 2ELsy 
_ - 2 
ee ee) Sey Jee! F547 (2) 
ae ae Giga | ea 
it+l al a lds Ue wal 2EIs ay OE. 4 KGA. 44 
3) 
ag Vee AG (4) 


It is to be noted that the mass term m includes the mass of 
hull and cargo and the entrained water or "added mass." 
This will be discussed in Section C. 

The relationshipSbetween the loading (shear), bending 
moment and deflection curves are well known from strength of 
materials. As w increases, the effect of including the rotary 
inertia and shear deflection terms in the equations becomes 
obvious, the shear and moment terms are altered increasingly. 


The boundary conditions for this free-free beam are 


OLS, 





and 


where subscripts s and B represent the stern and bow respec- 


tively (more accurately the after terminus and the forward 








terminus). Thus for station one 
Vi = Ve = 0 (5) 
ML v.v* 
2) = oP Se - ae me (6) 
S 1 1 S 
M, = M_ + Vi,L,.-J,w 9, = ~Jw* 6 Ce) 
1 S | ini Eee i S 
Z 5 
vo = spo ee wsl fe ‘oo ee 
il S sl 2EI, GEI, KGA, 
og 
va. = Vs ss 9 hy (8) 


It can be seen by inspection of equation (1) through 
(4) that the values for the bending moment and shear force at 
any point in the ship will be a function of an oe and w in 
addition to the calculated ship parameters. If, at the start, 
M is set equal to one, a unit defection of the stern, then 
the equationsfor the shear and bending moment at the bow will 


be expressed by the following: 
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V. = a + bO 
Ss 


M, =e) op de 

where a, b, c and d may be functions of w. From the boundary 
condition requirement it is known that both of these equations 
must equal zero. Arbitrarily setting M, = 0 and Ve =at 2s. 
and performing the sequence calculation a number of times for 
various values of w a curve of Ve versus wW can be obtained. 
This curve represents the exciting force Ve which would need 
to be applied at a specific frequency w to generate a unit 
deflection at the stern. The zeros of this curve, that is 
where the curve crosses the w-axis, represents the natural 
frequencies of the vessel where, with the absense of damping, 
a vanishingly small exciting force would set up synchronous 
motion. With the natural frequencies obtained the governing 
equations (1) through (4) will yield the displacement, shear, 
moment and slope for each station. The displacements for 
each frequency Ww) are considered for this discrete case to 

be equivalent to the continuous system eigenfunctions and 

are called eigenvectors. Indeed as the number of stations 
becomes large the eigenvector will converge to the eigenfunc- 
tion. As for the continuous case, the eigenvector must be 


normalized. A convenient normalization scheme is according 


to the following 
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N 

> M.y. = 1 ce) 
a 

which is the discrete representation of (32) in Section A. 


C. Added Mass 

The preceding technical treatment of the free and 
forced vibration of a ship like beam involved a mass term 
which was either uniform or not along the length. This term 
included not only the mass of ship structure, machinery, 
cargo, etcetra as might be obtained from a weight curve but 
also a hydrodynamic "added mass." When a vibrating body is 
immersed ina fluid (air or water, for example) the fluid 
Surrounding the body is put into continual motion which 
requires the expenditure of energy. Under conditions of 
synchronous or resonant motion, the body with this entrained 
water vibrates. This additional mass serves to lower signi- 
ficantly the natural frequencies of the vessel. It is 
therefore critical that a determination of this mass be made 
and added to that of the vessel. The total mass then, called 
"virtuai" mass is the sum of the ship mass and the added mass. 
This phenomenon has been observed and treated for a few 
Objects in classical hydrodynamics. In the case of an 
infinitely long circular cylinder having a cross-section of 
radius r translating perpendicularly to its axis with velocity 


Vin an infinite inviscid fluid it is possible to derive the 


Bz 





potential flow and to calculate the kinetic energy in the 


mPa. it is found to be [20} 


per unit length of cylinder where pop is the fluid density. 
Letting M' = om it is apparent that M' is the mass of 

mera that is displaced by unit length of the cylinder. It 
1s also apparent that the effect on the cylinder motion from 
the presence of the fluid can be allowed for by an addition 
of M' to the mass per unit length of the cylinder. An energy 


balance equation for forced motion becomes 


where F is the external force. 


OnesimplLiryving 
av 
' —_— = 
(M+M ae 1c 


or 
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Mere the presence of the fluid is equivalent to a force per 
unit length opposing motion. Obviously this virtual mass 
effect is only present in the case of accelerated motions 
such as vibration. For a sphere of radius r the kinetic 


energy of the fluid has been shown to be [20] 


Here M‘' = Sonr™ which is one half of the mass of displaced 


fluid. In general 
M' = C times the mass of displaced fluid 


where C is a vertical added mass coefficient. Comparison 

of the added mass coefficients for the cylinder and sphere 
both of radius r shows the former to be twice the latter. 
This reveals the three dimensional effect. The general 
method for determining the added mass coefficients of a 

Ship hull is a stepped procedure utilizing first a two 
dimensional analysis then multiplying by a three dimensional 
correction called a "J" factor which will be discussed later. 
The discrete analysis of the Prohl sequence is inherently 
compatible with this. At each station the geometric and 
material characteristics are assumed to be constant along 


the station length. Thus the station can be considered to 
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be a portion of an infinitely long cylinder and the added 
mass can be determined. The cylinder is supposed to consist 
of the underwater hull cross section and its above water 
mirror image. The actual ship added mass is thus one-half 
of this computed amount. This analysis is based on the fact 
that the fluid motion is streamline throughout, that is, 
irrotational. This is reasonable since the amplitudes and 
absolute velocities are small. 

This assumption of irrotational flow implies that the 
url of the fluid particle velocity vector field is zero, that 


1S, 


where V is the "del" operator. fThis implies that the 
velocity vector field is the gradient of a potential 


fomMction 


<I 
iI 


Vo (i) 


Additionally where the flow is incompressible and there are 
no sources or sinks a fluid continuity equation can be 


written in the form 
- aV 


Sys) 





memudoon Substitution of (1) into (2) 


Vo = 0 (3) 


This is the well known Laplace's Equation. The kinetic 


energy of any flow where a potential exists is [20] 


oT = f p(¥) “av 


where dV is a volume element and ep is the fluid density. As 
has been shown the added mass can be determined from the 
kinetic energy of the flow. It remains then to compute the 
potential for the flow around a ship Section. 

Classical hydrodynamics has evaluated the potential 
for a very few geometric shapes. In two dimensions flows 
past circles and ellipses have closed form equations and 
spheres and ellipsoids of revolution have equations for the 
potential in three dimensions. A method which greatly 
expands the number of shapes which can be analyzed is called 
“conformal mapping" which involves application of analytic 
function theory and functions of a complex variable. 


Tomer lustrate, let 


w= £(2) (4) 
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where 


see 


Za ee ay, 


(4) then represents a mapping of all points in the z plane 
onto the w plane. If the function f(z) is single valued 
then the mapping is one to one, that is to each point z 
where £(z) is defined there exists one and only one value 

of w= u+ iy in the w plane. The inverse will also be 

true if f(z) is analytic and f'(z) # 0. If these conditions 
are met, then it is also true that relative angle and shape 
are preserved in such a mapping (thus the term name "con- 
formal"). Right angles in the z plane will be mapped to 
right angles in the w plane. Another property of complex 
analytic functions is that the real and imaginary parts 
satisfy Laplace's equation and therefore can represent fluid 
flows. Knowing that the potential and stream functions 


Satisfy the Cauchy-Riemann equations 


oo ~ ov 
OX Jy 
Leo) 
oY OX 


OW 





and Laplace's equation (3) an analytic complex function can 
be generated by setting the potential equal to the real part 


and the stream function equal to the imaginary part 


(2) = $(x,y) + iv(x,y) 


This function is often called the complex potential. Now 


returning to (4) and writing its inverse 


Z = F(w) 


it iS apparent that the complex potential $(z) is expressible 
as a function of w in the form ¢[f(w) ] 

Further this function is analytic and therefore repre- 
sents a new fluid flow in the w plane. Each of the equipo- 
tential lines, boundary lines and streamlines in the original 
plane are mapped into the w plane. The orthogonal relationship 
is present and thus the actual flow around the transformed 
boundary is known. For example, the flow past a two dimensional 
body of arbitary shape can be determined (theoretically) by 
generating a transformation function w = f(z) which will map 
the points of the boundary of a unit circle in the z plane 
onto the boundary of the arbitary shape in the w plane. The 


inverse of the function is then substituted into the complex 
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potential for the flow in the z plane which is well known. 
Thus a complete solution for the flow about the arbitary 
shape is obtained. In practice the inverse may not be 
manageable or only a rough approximation to the desired 
shape can be transformed from the unit circle. cCross- 
sections of ships can, with good accuracy, be transformed 
from a unit circle however. 

The first and definitive work in this area was pre- 
sented in 1929 by Professor Frank M. Lewis in the Transactions 
of the Society of Naval Architects and Marine Engineers. The 
language of the paper is a bit quaint and some of the ques- 
tions raised have long since been answered, but the method 
presented for computation of added mass has seen only very 
minor modifications to the present day. 

Let x,y be the coordinates of any point in a two 
dimensional flow and ¢, the values for the velocity potential 


and stream function at the same point. The equation 


Z=x +t iy = £(o + ib) 


where £f designates an arbitrary function represents a two 
dimensional irrotational flow. The flow in the x, y plane 
can be further transformed to any other plane X, Y by the 


relation 


Z=xX+iye f(x + iy) 
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where ft. Tomanocierearolerhatry Eunction. 
Still a function of $ + iv. 


Take the flow past a unit circle 


and transform it by means of the relation 


Z = X+1Y = oe = eg ae 
Z 

Now 

Ze eh yy = rei? = any 
sO 

X + 1Y = me - nea ~ nee ee 
Since 

a = cos§ + ising 
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Z is of course 


soe 


(6 ) 


Cf 





it follows on equating real and imaginary parts that 


(l+a) cos8 + bcos398 


~ 
II 


(8) 


Y = (l-a) sin? - bsin38é 


These are the parametric equations of the transformed circle. 


The semi-axis parallel to the flow (8@=0) will be of length 


Y, =l-a+t+b 


femthne ratio of half beam over draft will be 


l-at+b 
lt+tat+b 





H = B/D = 


For the flow past a unit circular cylinder $¢ and yj are given 


by 
1 Ib : 
d6 = U(r + =) cosé yw = U(r - me Sing 


and on the boundary n=l 
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¢ = 2Ucosé y = 0 ee, 


This represents the flow at the boundary of a stationary 
cylinder and, also represents the flow at the boundary of 
the transformed cylinder. To obtain the solution for a 
meveng cylinder in a fluid stationary at infinity a super- 
position of a uniform flow of -U must be made. For this 


Mawrtorm flow 


@ = -UX and w = -UY (11) 


so uSing (11) and (8) in (10) 


U[{(l-a)cos98 - bcos36] (ae) 


-©- 
I 


-U[(l-a)sin9 - bsin36] (13) 


= 
I 


The kinetic energy of the fluid is 
2T = oss ive a. Vv") dxdy 
x y 


and since 


= B04 9) 2 
2T = psSL(GE)” + (58) “laxdy 
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An application of Green's theorem allows the double 
integral to be replaced by a line integral computed on a 
path which encloses the area. In this case the fluid is 
bounded by the cylinder surface and by a circle whose 


radius is allowed to approach infinity. 


27 = pM (ze)* + (ge) “ldxdy = - pf og% ds (14) 
ore aC 
a 2 


where ds is an element of path length, 3n is an inwardly 


pointing element normal to the curve, Cy 1s the cylinder 


circumference and C5 represents the circle with radius 


approaching infinity. Since the fluid is at rest far from 


the cylinder the line integral over Co will be zero. Since 
ST al) 
on ds 


the kinetic energy is given by 


2T = pfoap (15) 


Al 
which considering that the expressions for ¢ and w are 


dependent only on 9 equals 


8=2T 
2n-= pf U[ (l-a) cos8-bcos386] [-U[(1-a) cos§-3bcos36] ]dé 
0 


Tis Gia Oe (16) 
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mee the Semi-axis Of the transformed circle is to be B 


instead of (l-atb) this energy becomes 


Be 


mU“p[ (1-a)*#3b*] ——2 —, 
eae 


a 


2, (1a) *+3b* 


TUB 5 
(l-at+b) 


J er 


As previously mentioned the kinetic energy of flow past a 


erecular cylinder of radius B is 


therefore the virtual mass coefficient for the transformed 


@irccle is 


2 2 
a (l-a L+38 (18) 
(l-atb) 
and (17) becomes 
27 = = M'u" 
T = C7TB 0PU = MU 
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faeeene sham hull floating on the surface this is reduced 


by one-half so 
eee 
Mt = SCmB'p (19) 


Solving (18) and (9) simultaneously for a and b yields 


~ (2CH-H+1) +(H+1) /4C=3 


4 Sa 
2 (CH*-H*-H-1) 


] (20) 


Se ebaa ey, 
a= i (21) 


and the parametric equation of the transformed circles when 


of draft D instead of (ltatb) are 


x p pitta) cosétbcos 36 


= (l+atb) ee 
7 (l-a) sin8-bsin38 
5G 2 ar aera (23) 
X in this case being vertical and the line §@=0 coinciding 
with the X axis. [In [22] is a collection of graphs for 
Various values of H(half beam to draft ratio) and C. The 
H values correspond to common ship cross sections. For each 


of the paired values of C and H the constants a and b were 
computed from (21) and (20) then the shapes of the transformed 
Circle were plotted from (22) and (23) for values of 6 from 

0 to 7/2. The added mass for any ship section is then 
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computed as follows: a visual comparison of the shapes 
mri22) is made with the body plan for the station. Of 
course, the body plan must be properly dimensioned and the 
value of H known. The C value for the actual ship section 
is obtained by interpolation between adjacent curves and the 
added mass per foot for the ship section is obtained from 
(15) remembering that B is half beam. This procedure is the 
basis for present day calculations. Some refinements which 
make the procedure programmable are presented ina paper by 
L. Landweber and M. Macagno in 1967 in the Journal of Ship 
Research. An equation of the form 


A A 
(24) 


(of which (1) is a truncated version) can describe either 
mie mapping of a hull shape to a circle or a circle to a 
hull shape. The constants A. will of course be different. 
Furthermore, inversion formulas exist which enable the co- 
efficient of one transformation to be determined if the 
inverse coefficients are known. Let (24) represent the 
mapping of a wetted hull-shaped contour consisting of a 
wetted hull and itsS mirror image across the waterline from 
the z plane into the w plane. S is the bounded area of the 
moet contour and Xo is the radius of the circle in the w 


plane. An expression for the added mass coefficient C. aS. 21] 
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_ i ~~ 
Cc. = tty = (25) 


Let the inverse of (20) be of the form 


B B B 
eee 3 
3 

Z Z 


p= Za (26) 


= 
A method of determining the coefficients B. is based on the 
property that, among the closed curves obtained from (22) 

by varying the coefficients, the curve which will bound the 


Maximum area will be a circle. Let the ship curve be 


termed Cr and its mapping in the t plane be termed Ci: ilies 
"2 = r(g)er? then 
d Pecize—"ommin( Gir irde) = io r*de Agate 
Cp Cp Cp 
where z is the complex conjagate of z. Simlarly for Ci 
6 dt = Zis, (27) 
“Al 
where Si is the area bounded by Ci: Also 
B B B 
Tee ee es yg 8 
ee oa) a = =F —e ee (28) 
Z Z Z 
and 
B 3B 
= 1 ee 
abt = (1 5 ee ya Z (2 9)) 
Z Z 





Atter substitution of (28) and (29) into (27) the condition 


that Si be a maximum that is 


yields a set of n linear equations for the coefficients By, 

i =1ton. Standard matrix manipulation techniques are 
available for then determining the B,- The difficulties 

with this method are that the convergence rate with respect 
to n is relatively low and the accuracy of some intermediate 
terms decays with increaSing n. The solution to this problem 
is to cut off the first mapping at n=8 and then map the 
nearly circular curve Ci from the t plane into a much more 
nearly circular curve C. in the w plane where 

he rene 


tt 


The Cc. are determined in much the same way as were the B.- 
The coefficients A. of equation (24) can then be determined 
from the intermediate coefficients B. and C.- The added mass 
coefficient C., 1s then determined from (25). 

This thesis makes use of the above outlined procedure 
to compute the added mass for vertical vibration of hull 


sections as defined by offset data for each section. 
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Another significant item in this treatment is to 
calculate the three dimensional effect mentioned briefly 


earlier. Lewis, in his paper, talks of a "J" factor where 


Pe ccuc meme ete enemay Of Flurd 
Kinetic energy of fluid if motion 
is confined to transverse planes 


J obviously depends on the underwater configuration and will 
approach unity as the ratio of length to beam gets larger 
and larger. That is to say aS variations of the hull from 
section to section diminishes. 

An ellipsoid of revolution is one of the very few 
Shapes for which an exact solution has been obtained for the 
three dimensional flow. It is this shape that is used for 
the computation of the J factor for a vibrating ship. As 
the vibration mode goes up, and the number of nodes increases, 
it 1S apparent that the movement of water becomes more 
longitudinal and less transverse so that J. values, where i 
represents the jth mode, will decrease. The calculation of 
J. and J. was made in the paper by Lewis. The treatment is 
Similar to that presented for the two dimensional case in 
that an equation like (1) for the kinetic energy of the fluid 


1 used 
—_— oe ad 
Zo eS Sos ds (eae) 
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where in this case ds is an element of area. The solution 
of this equation can be obtained due to a closed form solu- 
tion to Laplace's equation having been found for this flow. 
That is, the potential ¢ for vibration of two and three nodes 
can be substituted into (31) and values for the actual three 
dimensional kinetic energy can be calculated. The denominator 
in the expression for J (30) is computed by summing (integra- 
ting) the values obtained from a two dimensional analysis of 
cross sections of the body. Lewis tabled the values of J. 
and J, for a variety of ratios of L/B. The values are 
Shown in Table l. 

The analysis performed by Lewis required that the 
Vibrating cross section remain in the same plane, that is, 
the motion results from pure shear. This does not reflect 
reality and consequently there has been much discussion on 
the subject beginning with a written comment of Lewis' paper 
by his peer J. Lockwood Taylor [22] [30] and continuing to 
the present day. 

Current procedure [32] involves the determination of 
J by accurate measurement of the natural vibration frequencies 
of a given body in air and in water. For Lewis' circular 
ellipsoid of revolution a simple analysis is possible because 
the distribution of added mass is directly proportional to 
the distribution of actual mass (Cc. = constant along the 


entire length). The mass distribution is 
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where r is the radius of the circular cross section and a 
memetlon Of longitudinal position, 9 is density and subscript 


B refers to the body. The corresponding added mass is 
M' = C gnr“p (33) 
V W 


where Cal and subscript w refers to the water in which the 
body is immersed. J here is considered to be a constant 
value applied at each cross section so that the actual three 
dimensional added mass is obtained from the two dimensional 


analysis performed. From basic vibration theory [22] 
aes VK/M 
where K represents stiffness and M mass. Put another way 


(34) 


where £ represents frequency and subscript a represents air. 


This equation for the circular ellipsoid considered is 


aul: 





LEWIS 
L/B Bos 
00 1.000 
30 
D2. 4 .954 
16 
ips). 82 .922 
12.95 897 
ez 
io. 02 .859 
9 
8.48 .829 
Foal .789 
7 
5.83 STi 
5.07 .704 
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TABLE 1 
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1.000 
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859 

837 
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2 zZ 
Fa,2 - P,/Tr dL+3p Jrr aie 
Be > i= 27, | 
fo Pair dL 
JO 
a (35) 
PB 


where the vibrating mass in air M. 1s assumed to be due only 
to the mass of the body and dL is differential length. 


mguation (35) 1s reduced to 


Ui (ey) ee | (36) 


Values from this equation obtained by R.L. Townsin [32] are 
listed in Table 1 and a comparison with the Lewis values 
shows a significant decrease reflecting the fact that no 
restrictions or assumptions were made concerning the cause 
of the motion. For ship-like vibrating bodies the procedure 
is somewhat more complex. Frequencies in water and air are 
both measured experimentally and calculated. For the water 
calculation several assumed values of J are used to obtain 
a plot of oF versus £. A reduction constant obtained from 
comparison of the two air frequencies is used to correct 
the measured water frequency. This abscissa is used to 
obtain the J, ordinate from the ae versus £ plot. Townsin 
proposes and this thesis employs the following empirical 


mermila for J as earumecelon Of B/L 
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J 4 ie — oe IN) 8 / i, eae. i (37) 


This formula agrees closely with data through Ja. 
The remaining consideration which is presented deals 
with a "sectional longitudinal inertia correction" labelled, 
by the author, Ws (the lower case to reflect a relation with 
the Lewis J factor). In most determinations of actual added 
mass per section an equation of the form of (19) is consis- 


tently used. That is 


ee 2 
m+ = 5C, 7B (38) 
and 
t 
Me = 
Dee 
1 
Or 
N 
M' = J ¢f m. (39) 
aL 
i=l 
so 
ear Z 
Mi = 5C;B.pd (40) 
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where Mm. represents the 2-D section added mass, C. and B. 

are the section added mass coefficient and half beam respect- 
ively, M' is the 3-D added mass for the entire body and Ms 

is the supposed 3-D section added mass. The departure from 
reality in (40) is apparent. J of the entire body represents 
a fixed fraction which when multiplied by the added mass of 

a 2-D flow gives the actual 3-D added mass for the entire body. 
There are some sections of a vibrating ship where the fluid 
motion is in fact very close to two dimensional and the added 
mass should be very close to that calculated from (38). At 
other sections, particularly near bow and stern, the flaw 
departs quite significantly from a pure two dimensional flow 
and the added mass would be less than that obtained from (40). 
To take into account these aspects the author proposes to 

use Jae 


Let the added mass at section i for a 2-D analysis be 
Ass 2 
= 5C, 1B. 0 (38) 
and let the actual 3-D added mass be 

Z 


fi» °: 
“gia ee oe 


5S 





then 
N N 
Mi = 5 M" = f ole aie (42) 
i - 3 


The requirement that equations (39) and (42) yield the same 


value for M' requires that 
N 
= Te (43) 
2 
Let 
jee = a = deb (44) 
a 


where d. is a difference factor computed for each section by 


N 
1 ee fs m2 mr 
Al. c= Ree Cae aa | eS 
i D oN 
N 
T 
[my 7Mz 1! ° = ei ite 
= 2 2 SO i43} 


and b is a normalization constant to be determined by substi- 
fmetom Of (44) into (43). In (45) N is the number of stations 
and n is the number of nodes in the vibrating mode. As the 
number of nodes increase three dimensional flow increases at 
all stations so J, 79. An example will serve to illustrate. 


Assume that an analysis has determined the 2-D sectional 
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added masses to be as listed in column 2 of Table 2. The 
ppckerence Eactors trom (45) are listed in column 3. Sub- 
Stituting the totals of columns 2 and 4 into (43) yield 


Hoist 2 eGo Ze 2o500 + 47886n) 5 


bese, (osc 47 6o6n) 


nie Oe) a 
4.633 0 
1.723 1 
1.059 , 
0.764 3 


Column 5 lists M. Gig  eoxs) mJ. Columns 6, 7 and 8 show the 
actual added mass at each station for vibrations of 0, 1 and 
3 nodes respectively. Notice as n increases the added mass 


approaches the values in column 5. 
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D. Wave Excitation 

Ties soreing Lunctlon F(x, ,t) used in the previous 
treatment represents the total external force per unit length 
at any section Xs and any time t. It includes the many effects 
of the waves encountered, the motion of the ship and the fluid 
flow. It is, as such, very complicated and the subject of 
much research. [12] [16] 

The most simple formulation of wave excitation would 


be due strictly to buoyancy changes from wave passage. 


where an elemental volume dV is equal to b.dy to give 


Ce) 
E(x, ,t) = ‘ yb, dy (1) 


where b. is the breadth of the ship and a function of y, 

Cc is the wave height at any time t and y is the specific 
weight of water. It has been traditional in naval architec-~ 
ture to use a trochoidal wave [14] to model the wave form 
for deep water waves. With sharper crests and flatter 
troughs it makes a very close approximation to the actual 
sea surface. However, this wave form cannot be derived from 
the velocity potential for surface waves and the parametric 
equations which describe a trochoid are an additional 
complexity in any analysis. 


Sys, 





A wave form which is more manageable and is derived 
from hydrodynamic theory of waves of finite amplitude was 
prepared by G.G. Stokes in 1847 [14]. 

Pe Kee 

c= Ecos (wt-Kx) +-5— cos2 (wt-Kx) (2) 
where t is half wave height from trough to crest and the 
Circular frequency w and the wave number K are related by 

K = oe = 2m/L, =" G/ Vv (3) 
for deep water waves where L is the wave length and ve is 
the celerity or wave velocity. In (2) a simple traveling 
coSine wave is modified by a harmonic which is half the 
length of the fundamental but which moves with the same 
celerity. This wave form has the same sharp peaks and wider 
troughs as does the trochoid but the fact that it is nonlinear 
with respect to height makes it somewhat less attractive. In 
the program suitable approximationSwill be made. 

It is well known that the pressure in a still body of 
water is equal to yh where h is the distance below the still 
Surface. It is also well known that in deep water the water 
particles have a circular motion which serves to reduce the 
pressure according to Bernoulli's theorem. That is, for a 


Simple harmonic wave, 
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[eee eens Tre “eos (wt-Kx) (an) 


This is the Smith effect. The pressure on the hull due to 
the wave passage is then 


p, =sye a? (5) 


where ¢ is given by (2). The buoyancy due to this pressure 
1S an integration of the vertical component of the pressure 
over the surface area of a hull slice of unit length. That 


is 


where n is an outward pointing unit normal to the hull and 
Mees the unit vector in the positive y direction. For a 


Ship vertically sided to the baseline this reduces to 


I 


18 (ena) p(T) b 


ies with (5) 


neds 


iD (epic) SS lovesaeam (6) 


where T 1s the draft. 
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Since the pressure given by (5) is seen to, in effect, 
reduce the waveheight f, an effective wave height ¢* can be 


defined so that 


ae 
¢ a 


£, (x,t) = ybr* (7) 


where the assumption is made that the ship 1s vertically 
Sided at the waterline so that b is constant over the range 
Seec*. The subscript 1 indicates that this buoyancy force 


Momonly part of the total excitation. ¢* is given by [16] 
me K : 
Be = G(l-- f ze “ay! (8) 
=-T 


where z is the hull offset and y is distance to the design 
waterline. The integral serves to weigh the determination 
of the reduction factor by the hull shape. For a section 


with a vertical side to the baseline (8) would reduce to 


0 CX 
c* = c[1-K f e ay] 


Al 
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so that 


£, (x,t) = yhee ~~ 


mime 1S 1dentical €6 (6). The discrete form of (8) is 


and for the discrete forcing function per unit length 


] 
= **— 2 aaa Y 


where f is given by (2). 


(10) 


The behavior of the bracketed term is to reduce the 


wave height as K increases. For small K the term is very 


close to 1. Consider the straight sided vessel again. 


reduction in wave height to one-tenth of the long wave value 


requires 
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For a draft of 20 feet K = 8.70 which corresponds to a wave 
Mength Of ./2 feet. So as K increases the correction term 
approaches zero. Thus there are two major effects taking 
place to reduce the energy available for excitation of the 
higher modes of a ship hull. First the wave height decreases 
as the frequency increases and the effect on an immersed body 
decreases as the frequency increases. 

Due to the logarithmic decrease in (9) the discrete 
form of the correction term, hereinafter called D, causes 
some convergence problems, that is,as K increases D fails to 
approach zero if the increment Ay. is too large. In the pro- 
gram four foot waterlines are used and the following slightly 


modified expressions must be employed. 


Wess: 
=e 


(11) 





Equation (ll) has the proper limiting values and 
accounts for the hull shape also. 

There are two other parts to the forcing function to 
be presented. As the wave passes down the side of the ship 
the vertical component of the orbital velocity has a damping 


effect and this force can be expressed by 


£ (x,t) = (N'-Vdm'/dx) 7* (2) 
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where (N'-Vdm'/dx) is the sectional hydrodynamic damping 
coefficient which will be discussed later and ¢* is the 
time derivative of the effective wave height. The third 


force is due to the inertia of the added mass 


£3 %57t) = fhe C* (23) 
where t is the effective vertical acceleration of the wave 
height. As was discussed in Section C the added mass may be 
incorporated into the equation of motion in two ways, either 
as an addition to the mass of the vessel and placed on the 
left hand side or as an exciting term and placed on the right 
hand side. The superposition principle allows the separation 
of the added mass effect in this case since two separate 
motions are being superposed, the vibratory motion of the 
vessel and the harmonic motion of the waves. 

ene ACS ) ae is the added mass calculated for immersion 
to the still water waterline. AS the wave passes, the draft 
at a section changes and so does the added mass, but here 
the sides at the waterline are assumed to be nearly vertical 
and the change in added mass is very small so the effect can 
be neglected. 


The total excitation force is then 


3 
dete)! £. (x,t) (14) 
cE 
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and the modal excitation force is obtained by multiplying 


by Y_ and summing over the ship length as shown in Section A. 
N 
F(t) = 2 ¥ (x, ) £ (x, ,t) (5) 
i=l 
The correct surface height, velocity and acceleration are 
_ DKE- 

ce Docos (wt-Kx) +=5— Co7 (or kx) 

c* = =Duwtsin (wt-Kx) -DKwr- Sin2 (wt-Kx) (16) 

a 2— 2-2 

6-5=— =p 2eCOS(wt-Kx)—-2ZDKy Cc Ccos2 (wt-Kx) 
and 

£(x,,t) =yb,c*+(N'-Vdm' /dx) ¢* + mo” 


= (yb, -m!D,w°Z) cos (wt-Kx) 


-D. wo (N'-Vdm' /dx) sin (wt-Kx) 


—2 
ees aay, 
+ ( 5 - 2D. Kw a m,) coS2 (wt-Kx) 
~D, Kwo* (N'-Vam' /dx) sin2 (wt-Kx) (17) 


66 





ret 


Gia = yb, Dc = m.D.wt 
Coo -D, wz (N'-Vdm' /dx) 
Sr. es. - 2m,KDjw 
Giy = -D, Kwo* (N'-Vdm' /ax) 
then 
£(x.,t) = G, 1, coswtcoskx + G,,Sinwtsinkx 


+ G...Sinwtcoskx - G. coswtsinkx 


a2 Zz 


-- G, ,cos2wtcos2Kkx + G, ,Sin2wTsin2Kx 


cw SinzwecosZkx = Gs,coszZwTSinZ2Kx 
14 14 


e@eseitution into (15) yields 


N 
F (t) = coswt > WOES aay) 2G. 
r aL ak 


cosKx A 
a i=l 


il 


N 
+ Sinwt 2% VowixG..sinkx.B 
a es er = ile 


oe 
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Sine yoetx.) G..CcOSKxX.C 
eA: eZ aL 


coswt 2 a) G, 5Ssinkx,D 


Coszuc 2 Y(x.) G..cOSZKX. 
me J, 18: al, 


Simczwe © YY (xX.) G..sin2K~x. 
iy iS a 


Sic, Yo 1s.) G. ,cOSZKx. 
ee 14 aL 


GOsvime. Y (X.) 46.,51n2Kx: 
emis al a 
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11 COS (wt-8,) 


r2C0S8 (wt-S5) 


races (2wt-8,) 
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where 


N N 
a 2 2/2 
Pia es Y¥. (X;)G;,Ccoskx,) cm Y. (x; )G,,Sinkx, ) } 
i=1 i=1 
N N 
_ 2 ; Bay 
P 9 = oe Ms (x, )G, ,coskx, ) +(Z Y. (x,)6; ,sinkx; ) } 
i=l i=l 
(22) 
N N 
Pee RE) Chee oar.) f(D v(x. )G. .sindkx.)-}/- 
FS ee | lc qaie: ais Re: vt Cpe tas ao on ear a att 
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N N 
i 2 ps, Maly A 
eA {(2) Y, (x;)G, ,cos2Kx, ) a4 Yr (x;)G,,Sin2Kx, ) } 
1=1 i=l 
N 
is eae nae 
Oo = arctan xq 
rr a aE a 
N 
Pies secre gees 
O45 = 7r-arctan y 
Hecate 3 
(23) 
N 
Pee io Ge ata 
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The ay i-l1,2,3,4 are known as "participation factors" 
and were obtained by weighting the wave forces by the mode 
shape. The factors indicate the participation of the wave 
forces in exciting a vertain mode. Solution of the independent 


equationsfor the modal response as discussed in Section C 


yields 
4 P.,cos(wt-9,+¢e_) 
n(t) = 2 (24) 
7 i=l é Ray ome eee aa 
sag 1g 

where ee 

- = arctan ao. 

=O, 


The most significant part of this exact time history 
is the magnitude of the response. Equation (24) may be mani- 


pulated to yield the amplitude of the sinusoidal response 


4 4 
_ 2 Ph lyf 
fla = {(2 P.,cos8.)" + (Z£ P ,sind,) } (25) 
i=l i=l 


The nonlinearity of the expression for the height of 
the sea surface can be eliminated by dropping the modifying 
harmonic. Although the simple cosine term is not an exact 
representation of the sea surface it is ideal for representing 
the many component waves of the general sea. It also allows 
the principle of superposition to be applied and for these 
reasons the program encorporates this modification. Result- 
antly G,=G,=P,=P,=0 with significant subsequent simplification 


5 oe ee 


of the program. 
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This analysis assumes that the presence of the vessel 
Siees not Gisturb the wave motion which it, in fact, does. 
This, then, provides an upper bound to excitation and real 
forces will be less. How much less is a matter for speculation 


and research. 


E. Damping 

Tt would not be incorrect to say that relatively little 
is known about damping associated with ship hull vibration. 
To be able to calculate the amplitudes of desired quantities 
such as stress levels, bending moments and displacements, 
the distribution of damping in the hull, cargo and surrounding 
water should be known. The available literature is sparse, 
[4] [11] [16] [17] but it is generally assumed that energy is 
dissipated thru a number of mechanisms. A convenient cate- 
gorization divides the mechanisms into those which are hydro- 
dynamic and those which are not. Water friction and the 
generation of both surface and pressure (sound) waves are 
the hydrodynamic effects. Structural and cargo damping are 
non-hydrodynamic. Another possible way of differentiating 
the loss mechanisms is to determine where the damping takes 
place, either within or external to the hull. 

As described previously, superposition of the natural 


modes of the ship yields for vertical vibration 


viet) = EY (x) my (t) (1) 
tal 





where Y (x) is the normalized mode shape and n(t) is the 
natural coordinate which incorporates the natural frequency, 
initial conditions and modal damping factor into a time de- 
Pemdent function. Motion for r=0,1 1s heave and pitch 
respectively. Fortunately the damping in ship hull vibration 
is relatively small so that calculations for the natural 
modes and frequencies can be carried out assuming no damping 
without significant error. This is especially true for the 
lower modes. 

Other qualitative information concerning damping has 
been known for some time. Higher modes are less affected by 
hydrodynamic damping due to in part the smaller displacements 
involved. In these modes structural damping is the major 
factor. The situation is reversed for lower modes. The 
generation of sound waves increases aS the vibration fre- 
quency goes up but this remains a very small quantity through- 
Out the range of interest, that is w less than about 35 rad/ 
Sec. For rigid body motion there is, of course, no structural 
damping. If vibration analysis is of concern and sea-keeping 
is not then discussion can be limited to distortion modes 
ae 2. Surface wave generation increases with decreaSing Ww. 


For even the largest of vessels, however, ls greater than 


ay. 
about 2.5 rad/sec and wave making damping is far less than 


Sebuctural damping. [17] 
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imc Giscussion Of Gamping Lt 1S common to refer to 
the logarithmic decrement 6 associated with a mode or given 


damping mechanism. 
oo 2nc/Vv1-c2 = 1n(x,/X,) (2) 


where t is the critical damping ratio associated with the mode 
or mechanism and xy and XxX». are values of the response separated 
by one wave period. 

Since structural damping is the significant factor and 
appears to be manageable there has been some effort directed 
at its analysis.[{4] This damping derives from hysteresis in 
the steelwork. The logarithmic decrement for this is a power 


function of the stress. For levels of stress less than about 


30,000 psi 


where C is a material constant and n = 3. At higher stress 
the exponent increases rapidly and other effects such as 


stress history become important. For steel plate (3) becomes 


1) 





Values for a welded structure go up with the complexity of 


the structure and may lie in the range 


This dramatic increase is explained in part by the high 
residual stresses from welding and stress concentrations at 
joints. There may also be as yet undiscovered hysteresis 
effects which are structure related instead of material 
related. 

Ships of current design are usually welded but in 
riveted ships coulomb friction at the joints is very important. 
In traditional analysis this dry friction is treated differ- 
ently from other damping effects. The coulomb damping force 
is not a function of velocity rather it has a constant value 
#Or all velocities and varies only in sign. This directional 
force is treated as an excitation force which opposes motion 
and equations for each direction arc solved. One test 
concluded that riveted ships have twice the structural 
damping of welded ships.[17] This would obviously vary 
with ship age and history. 

Experiments have been performed on loaded and unloaded 
vessels to determine cargo damping effects. Results show 
that cargo does increase the logarithmic decrement but for- 
mulations which have been made to quantify the effect have 


failed to predict well. [17] 
74 





The literature separates the generation of surface 
waves from the other damping mechanisms. It is that part 
of the hydrodynamic damping that can be derived from potential 
theory and the only part of ship damping with a solid mathe- 
matical base. However, the literature goes on to say that in 
the distortion modes the natural frequencies are relatively 
high and the displacements so small that for most existing 
vessels the effect can be taken as zero. A bit ironic. In 
equation (11) of Section D the term in parentheses (N'-Vdml/dx) 
is the sectional hydrodynamic damping coefficient according to 
the strip theory approach used in this paper. N' is a para- 
meter which quantifies the generation of surface waves and 
Vdm/dx is a forward speed correction term. V is velocity and 
m' is the added mass. This coefficient, it 1S seen, varies 
along the hull length and this distribution of damping effect 
is important for computation of exciting forces. The normal 
mode analysis in Section A required that the damping be pro- 
portional to the mass and stiffness distribution. Otherwise 
the modal equations would not be uncoupled. Since the pre- 
dominant damping SOuUpeeCeroestrUCEUral this aSsSumption 1S not 
too far off, even though the distribution of the sectional 
hydrodynamic damping 1s weighted at the bow and stern. N' is 
nearly zero amidships but gets larger at the bow and stern. 
The effect is still very small. Only the forward speed 
correction will be considered in this work and consequently 
the amplitude obtained will be slightly high. 


es 





Gemeotaehattonwer equation (2) brings up the fact that 
in much practical usage all the damping effects are lumped 
into one parameter and its value obtained from a spectral 
analysis of stress data taken from a vessel which is loaded 
impulsively so that all modes are excited. Several full scale 
measurements from slamming or mechanical excitation (anchor 
drop test) have been obtained and modal damping constants 
calculated. [{11][17] Some convenient empirical relations have 
been devised from such experiments and other means so that 
values for the C.. can be obtained quickly. From basic vibra- 


tion theory [24] 


/ (21) *+82 


and a formula for the logarithmic decrement which seems to 
give intermediate and hopefully representative values is 
a7 J 


$§ = 1.065 x 10. GE'S Gay, 


16 i oe 


Equations (4) and (5) will be used in this paper to determine 


mae C . 
1g 
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ae PROGRAM OVERVIEW 

Estimation of the parameters for use in the Prohl 
sequence calculations is a choice made by the author. The 
values of m, A, I and J for each section could be obtained 
from scantling lists and a weight curve but their assembly 
would be time consuming and defeat one of the purposes of 
this work, that is rapid vibration information with a minimum 
of data required as input. Therefore the subroutine HULL for 
estimating these quantities has been incorporated into the 
program. Input required is station spacing, number of 
stations, displacement, thickness of side shell, main deck 
and bottom plating, offset data and deck edge heights for 
each station. 

Sectional mass is computed by first obtaining the 
total sectional area from the offset data for each station. 
A partial sectional area is computed by the use of the five- 
@€ilght rule [14] and the offset data. An additional area is 
computed which is equal to the deck edge height minus the 
highest waterline for which offset data 1S input times the 
offset at that waterline. The sectional area equals the 
sum. A uniform permeability is 4&4 reasonable assumption for 
naval vessels so that the mass is proportional to the area 
at each station. The sum of the sectional areas is divided 
into the displacement to get an area density. The density 


is then multiplied by each sectional area to obtain a 
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sectional mass. No attempt is made to account for super- 
structure concentration of mass or its effect on stiffness. 
A two beam ship model would be the next step in sophistication. 
The mass moment of inertia, J, about the neutral axis 
is determined by first assuming that the beam cross section 
has a uniform distribution of mass. Naval vessels are 
characteristicly longitudinally framed. In the midship 
section the hull strength is provided by the main deck, 
Side shell and bottom plating and their associated longi- 
tudinals. A thirty inch spacing is not uncommon and the 
stiffener area to plate area ratio iS about .3. So, the 
effect of longitudinal is taken into account by increasing 
the plating thickness by 30%. The neutral axis is determined 
by a calculation of the first moment of the hull steel. The 
sectional area density is then used along with the offset 
data in a second moment calculation about the neutral axis. 
Each hull cross section is divided by waterlines with four 
Foot spacing. The area from the baseline to the first water- 
line is assumed to be bounded by a parabola. All subsequent 
areas are assumed to be trapezoids. The distance from the 
neutral axis to the centroid of each incremental area is 
computed, squared and multiplied by the area. The second 
moment of the area about its own centroid 1s computed. The 
sum of all these terms equals the sectional second moment. 


This is multiplied by the sectional area density to obtain J. 


i] 
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The moment of inertia, I, iS obtained by a second 
moment calculation of the hull steel around the neutral 
axis. 

The effective shear area KA ina hull shape is very 
difficult to determine. Traditional practice has been to 
use only the vertical side shell and any continuous longi- 
tudinal bulkheads in the calculation. Analytically obtained 
values for K are for simple symmetric shapes only. [15] 
Calculations by the author for box type sections with 
various side, deck and bottom thickness and beam/draft ratios 
shows that an approximate value for the effective shear area 


can be obtained from the following 


(KA), = .87(2tD) = 1.74tD (1) 


where t is the sideshell thickness and D is the depth. Then 


mores tee | 


KAG . ene Dir 2 Gl Fi) 


I 


Or 


KAG; = 22S x. ae GE, (2) 


meet in inches and D in feet. 
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Subroutine NMODES first computes the virtual mass 
of each station which varies with the mode of vibration for 
which a solution is sought. Burrill’s formula {14] is used 
to estimate the frequency of two aeae vibrations and the J 
factor for that mode is formed and virtual masses computed. 
Equations (1) through (4) of Section B are solved for each 
station given a value of w in the following way: the shear, 
bending moment, displacement and slope at any station i are 


assumed to be of the form 


Meee en 4 Cd. G 
of a ivs 
> 
a ei Ee 
mea 4s 
where ore is an unknown initial slope at the stern. Then (1) 


through (4) of Section B are used to determine the conditions 


ae Station itl 
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From the treatment in Section B 
ee N 
S do (3) 
and 
lowe 
ay + = 0 (6) 
N 
at a natural frequency. The subroutine performs the sequence 


calculation until opposite signs are obtained for the quantity 
on the left hand side of (6) for two values of w. A natural 


frequency is then known to lie between these two values. 
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The subroutine is designed to determine the natural fre- 
quencies to an accuracy of +.04 radian/sec. After a natural 


frequency is obtained the mode shape is determined from 


y.=e.+f. = (7) 


Data obtained from the Andersson reference [2] for a 
60,000 ton tanker was input into NMODES. The resulting 
natural frequencies obtained are compared with the measured 


frequencies in radians/sec. 


NMODES (20 a 
1 4.92 4.89 + .6 
2 10.63 De 4s +1.9 
3 16.68 16.41 +1.6 
4 23.20 Dore Sl ee 
5 IN SS PSS) ANA See: 
6 34.88 B5..05 = ad 
7 Bot 40.33 -1.4 


The percent difference is quite small. 
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The input subprogram reads data which is punched 
according to the format listed after each read statement. 
The data must be placed in the proper fields, decimal points 
must be punched and all integers must be right justified. 
STA is an integer variable and is the number of stations 
for which offset data is given. STASPA is the station 
Spacing in feet, DISPL is the displacement of the vessel in 
long tons. The program later changes this value and does 
all calculations in slugs, BEEM is the beam of the vessel 
in feet. (Another variable with the correct spelling is used 
in the program.) Tl, T2 and T3 are the main deck, side shell 
and bottom plating thicknesses respectively of the midship 
section in inches. DECKE is the deck edge height from bow 
to stern in feet. There must be one entry per station. DATA 
is offset data, one card per station. All entries are in 
decimal feet. The first field is a real integer (decimal 
point must be included) which is closest to the draft at 
that station divided by 4. The remaining twelve fields are 
offset data starting from the bottom and proceeding upward 
at four foot intervals. The first field must have a real 
integer less than or equal to twelve and only 12 waterlines 


may be input, unused fields may be punched with zeros. 


83 





The output will list the computed station mass in 
slugs, the mass moment of inertia J in slug ft“, the second 
moment of area I in oe shear stiffness KAG in l1bf and 
the 2-D added mass for each station in Slugs. The program 
will list the first seven natural frequencies and the 
normalized made shapes. The program will also list the 
response amplitude operator over the frequency range of 
interest for speeds of 0, 10, 20 and 30 knots. It will also 
list the response of the vessel to a Pierson-MOskowitz 
fully developed sea with a wind speed of 50 knots. Response 
amplitude is in inches measured at the bow. 

It is obvious from a look at the Pierson-Moskowitz 
formula that precious little BORO erist s at frequencies 
above the first mode of vibration, springing. It is reason- 
able then to assume that wave excitation is significant only 
for large, long and shallow vessels. The large masses and 
the low stiffnesses (due to reduced depth) are responsible 
for a relatively low first natural frequency, one low enough 
Eadat SULLIcCient energy is available for excitation. Too, as 
wave frequency goes up the effect on a vessel of constant 
draft drops off logarithmicly as has been shown. The only 
times then that the higher modes will be excited are when 
Slamming occurs or under some other impulse type loading. 


The value of the program lies in the beam parameters obtained 
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from the HULL subroutine, the natural frequencies and mode 
shapes obtained from NMODES and the forcing data which can 
be retrieved from FORCE. The node location in the natural 
modes is important as possible locations for placement of 
machinery operating at or near that natural frequency since 
excitation at a node will not excite vibration. 

Numerous comment cards have been inserted in the 
Peecram and Can be used to understand the many small pro- 
cedures, which are a natural part of any project of this 
dimension. The text of this thesis serves only to discuss 
the major technical areas and not the interesting but 
Obviously minor structural details of software construction. 
The program listing which follows is divided into the sub- 


routines. A sample input record and output record are shown. 
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